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The relationships among charge order, spin fluctuations, and superconductivity in underdoped 
cuprates remain controversial. We use neutron scattering techniques to study these phenomena in 
Lai, 93 Sro.o 7 Cu 04 , a superconductor with a transition temperature of Tc = 20 K. At T ^ Tc, we find 
incommensurate spin fluctuations with a quasielastic energy spectrum and no sign of a gap within the 
energy range from 0.2 to 15 meV. A weak elastic magnetic component grows below ~ 10 K, consistent 
with results from local probes. Regarding the atomic lattice, we have discovered unexpectedly strong 
fluctuations of the CuOe octahedra about Cu-O bonds, which are associated with inequivalent O 
sites within the Cu02 planes. Furthermore, we observed a weak elastic (330) superlattice peak that 
implies a reduced lattice symmetry. The presence of inequivalent O sites rationalizes various pieces of 
evidence for charge stripe order in underdoped La 2 - 2 !Sra;Cu 04 . The coexistence of superconductivity 
with quasi-static spin-stripe order suggests the presence of intertwined orders; however, the rotation 
of the stripe orientation away from the Cu-O bonds might be connected with evidence for a finite 
gap at the nodal points of the superconducting gap function. 


I. INTRODUCTION 

The degree to which charge, spin, and superconducting 
orders compete or coexist is a topic of considerable cur¬ 
rent interest [1, 2]. In cuprates with spatially-uniform 
d-wave superconducting order, one generally finds the 
absence of spin order plus a gap in the antiferromag¬ 
netic spin fluctuations, comparable in energy to the su¬ 
perconducting gap [3-7]. In contrast, there are certain 
cuprates, such as La2-a;Baa;Cu04(LBC0), that exhibit 
some degree of spin and charge stripe order, where the 
coexistence of gapless antiferromagnetic spin fluctuations 
and superconductivity [8, 9] have been interpreted as evi¬ 
dence for a spatially-modulated pair-density-wave super¬ 
conducting order [2, 10, 11]. 

In the present paper, we examine two aspects of coex¬ 
isting antiferromagnetic and superconducting orders in 
La2-a;Sr2,Cu04 (LSCO) with x = 0.07 and supercon¬ 
ducting transition Tc = 20 K. First, we demonstrate the 
development of static spin order at temperatures below 
Tc, as well as the absence of any gap in the low-energy 
spin excitations. Study of the elastic scattering yields 
a detailed characterization of the structure of the short- 
range spin correlations, significantly extending previous 
observations on a crystal with x = 0.07 (Tc = 17 K) [12] 
and allowing an instructive comparison with x = 0.10 
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(Tc = 29 K) [13]. The spin order appears as incommensu¬ 
rate peaks about the antiferromagnetic wave vector, with 
the incommensurate peaks rotated away from the Cu-O 
bond directions, as observed previously in LSCO with 
X = 0.12 [14] and La2Cu04+y [15] (but different from 
LBCO [16]). We note that this is compatible with the 
identification of charge-stripe order developing at 90 K, 
as determined by a nuclear quadrupole resonance study 
by Hunt et al. [17]. The spin fluctuations have been mea¬ 
sured well below Tc and with several different incident 
neutron energies, providing reliable results between 0.2 
and 15 meV. The absence of a spin gap confirms and 
extends the results of previous studies on crystals with 
X = 0.07 [18], 0.085 [19], and 0.10 [20]. 

We also address the nature of the stripe pinning in 
LSCO. The atomic displacement patterns within the 
Cu02 planes for the relevant crystal structures are in¬ 
dicated schematically in Fig. 1. In the LTO phase, al¬ 
though the in-plane oxygens break C/ symmetry, their 
positions are nevertheless symmetry related, with all be¬ 
longing to the same Wyckoff position. In LBCO, it is 
the presence of two inequivalent O sites within the Cu02 
planes of the LTT (or LTLO) phase [21] that is associated 
with stripe pinning [22]; in both LTT and LTLO, the in¬ 
plane oxygens require two different Wyckoff positions. It 
is the resulting electronic anisotropy that can pin charge 
stripes which tend to orient along Cu-O bonds. 

Recently, the CuOe octahedral tilt fluctuations directly 
associated with the inequivalent O sites have been de¬ 
tected in LBCO by inelastic neutron scattering [23]. We 
have now detected these fluctuations in our LSCO sam- 
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FIG. 1. (color online) Diagrams indicating ont-of-plane dis¬ 
placements of the in-plane oxygens in the low-temperature 
orthorhombic (LTO), low-temperature less-orthorhombic 
(LTLO), and low-temperature tetragonal (LTT) phases; space 
groups listed under acronyms. -|- (—) indicates displacement 
above (below) plane; line thickness reflects magnitude of dis¬ 
placement. There are two inequivalent in-plane O positions 
in LTLO and LTT, but all are equivalent in LTO. 

pie, as well. In addition, there is a weak (330) superlat¬ 
tice peak that decreases gradually on warming to 300 K; 
its presence implies a lowering of the symmetry, likely to 
the LTLO phase, and a static inequivalence of planar O 
sites. The static symmetry reduction provides an expla¬ 
nation for previous indications of stripe pinning at rela¬ 
tively high temperatures in LSCO [17]. In addition, the 
relative intensity of the tilt fluctuations is large compared 
to the superlattice intensity even at low temperature, in¬ 
dicating the existence of significant entropy associated 
with these tilts. 

The rest of the paper is organized as follows. The fol¬ 
lowing section describes the experimental methods. The 
experimental results, including the characterization of 
the superconducting transition of the sample and the 
neutron scattering measurements of spin correlations and 
soft phonons, are presented in Sec. III. We discuss 
these results, their relevance to the concepts of spatially- 
modulated superconductivity and stripe pinning, and 
their connection to other recent work in Sec. IV. The 
paper is summarized in Sec. V. 


II. EXPERIMENTAL METHODS 

The crystals of Lai, 93 Sro,o 7 Cu 04 were grown at 
Brookhaven in an infrared image furnace by the 
traveling-solvent floating-zone method [24, 25]. After 
growth, the diameter of the cylindrical crystal was ap¬ 
proximately 6 mm. The superconducting transition 
was determined from a measurement of the tempera¬ 
ture dependence of the magnetization using a commercial 
SQUID (superconducting quantum interference device) 
magnetometer on a small piece of the crystal. 

For neutron scattering, two large crystals were cut 
from the original rod, with a total mass of 19 g. Prior 
to each measurement, the two crystals were co-aligned 
at room temperature to yield the total effective mosaic 
^ 0.8°, consistent with the « 0.5° mosaic of each of 
the two pieces. (For a photo of the co-mounted crys¬ 
tals, see inset of Fig. 3.) The domain structure of 


the low-temperature phase was characterized on triple¬ 
axis spectrometer HBl at the High Flux Isotope Reac¬ 
tor (HFIR), Oak Ridge National Laboratory (ORNL). 
The measurement was performed with a fixed final en¬ 
ergy of 13.5 meV and collimations of 15'-20'-20'-240'. 
Further measurements were carried out on the SPINS 
cold-neutron triple-axis spectrometer at the NIST Cen¬ 
ter for Neutron Research (NCNR). With the c axis verti¬ 
cal, scattering wave vectors Q = {h, k, 0) were accessible 
in the horizontal scattering plane. Wave vectors will be 
expressed in units of (27r/a, 27r/6, 27r/c) with the mea¬ 
sured low-temperature lattice parameters a = 5.324 A, 
h = 5.385 A, and c = 13.2 A. Based on neutron pow¬ 
der diffraction measurements, the crystal structure pre¬ 
viously has been identified as LTO phase [26]. 

Next, the sample was transferred to the Spallation 
Neutron Source (SNS), ORNL, where it was studied on 
the HYSPEC instrument (beam line 14B) [27]. In the 
transfer, a 2.7-g piece of crystal came loose and was re¬ 
moved, leaving a sample mass of 16.3 g. The sample 
was mounted in a Displex closed-cycle cryostat. Mea¬ 
surements were done with a fixed incident energy of ei¬ 
ther 8 or 27 meV and a chopper frequency of 300 Hz. 
The graphite-crystal array in the incident beam was set 
for vertical focusing. The energy resolution (full-width 
half maximum) near elastic scattering is ^ 0.25 meV for 
Ei = 8 meV and ^ 1 meV for Ei = 27 meV. For a 
typical measurement, the detector tank was placed at a 
mean scattering angle of either 33° or 80°, and then mea¬ 
surements were collected for a series of sample orienta¬ 
tions, involving rotations about the vertical axis in steps 
of 1°. (The position-sensitive detector collects neutrons 
for scattering angles of ±30° in the horizontal and ±7.5° 
in the vertical direction.) From such a set of scans, a 
four-dimensional (4D) data set was created and analyzed 
with the MANTID [28] and HORACE [29] software pack¬ 
ages. Slices of data corresponding to particular planes in 
energy and wave-vector space can then be plotted from 
the larger data set. For reference, the beam power was 
~ 0.85 MW, and the counting time for a single 4D data 
set was roughly 6 h. 


III. RESULTS 
A. Sample characterization 

The bulk magnetic susceptibility measured at low field 
on a piece of the LSCO x = 0.07 sample is shown in 
Fig. 2(a). As one can see, the bulk transition is at 20 K. 
The data have been corrected for the demagnetizing fac¬ 
tor, and within the 10%-uncertainty of that correction, 
the superconducting shielding fraction is consistent with 
100%. Measurements on pieces taken from different po¬ 
sitions along the crystal rod indicate a variation of Tj, of 
< 1 K. 

There is also some weak anisotropic diamagnetism 
above the bulk transition. To characterize that, the sus- 
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FIG. 2. (color online) (a) Bulk susceptibility of the LSCO 
X = 0.07 sample measured with a 20 G magnetic field applied 
along the c axis both for field cooling (FC) and zero-field cool¬ 
ing (ZFC) conditions, (b) Gomparison of magnetic suscepti¬ 
bility measured in a 1-T field applied either along an in-plane 
Gu-0 direction (open circles, labelled H\\a) or along the c axis 
(filled circles) after scaling by the anisotropic g factors taken 
from [30]. Inset shows the difference between these scaled sus¬ 
ceptibilities, indicating that diamagnetism within the Gn 02 
planes sets in below ~ 80 K. (To convert the susceptibility to 
SI units of Hm^ /mol, multiply emu/mol by (dvr)^ x 10“^®.) 


ceptibility was measured with a field of 1 T applied along 
the c axis and along an in-plane Cu-0 bond direction. 
Assuming isotropic spin susceptibility at high tempera¬ 
ture (> 150 K) and using the anisotropic g factors deter¬ 
mined in a previous study of LBCO [30], we obtain the 
corrected spin susceptibilities shown in Fig. 2(b). From 
the difference curve plotted in the inset, one can see that 
the diamagnetism of the Cu 02 planes becomes significant 
below ~ 80 K. This is consistent with previous results on 
LBCO, especially with x = 0.095 [31], and with studies 
of the onset of superconducting fluctuations in LSCO by 
torque magnetometry [32] and Nernst effect [33]. 

The twin-domain structure of the orthorhombic phase 
at T = 4.2(1) K is presented in Fig. 3, which shows mesh 
scans covering the positions of the (200) and (020) Bragg 
peaks. As indicated schematically by the arrows, we have 
two main twin domains. Least squares fitting indicates 
that the intensity ratio of the (200) domain to the (020)' 
is 1.08(6). 


B. Magnetic scattering 

1. Elastic 

We begin with the characterization of the elastic mag¬ 
netic scattering about Qaf = (1,0,0). A mesh of scans 
was measured at SPINS with Ef = 5 meV for a tempera¬ 
ture of r = 1.5(1) K, and the low Q part of the mesh was 



FIG. 3. (color online) Arrows indicate schematically the ori¬ 
entations of (200) and (020) fundamental Bragg peaks of the 
two main twin domains. Image plots show the results of mesh 
scans performed around these Bragg peaks. Upper right inset 
is a photo of the co-mounted crystals, each wrapped in A1 foil. 


also obtained at 30 K. As no magnetic signal was detected 
at 30 K, that data set provided a useful measurement of 
the background. The 1.5-K data were fitted with four 
gaussian peaks from each twin domain, at QApiqJ with 
Qaf = (1,0,0) and at ±91^ with = (0,1,0)', 
plus a gaussian in h and k to describe the background; the 
same background parameters were also fit simultaneously 
to the 30-K data. The 1.5-K data, after background sub¬ 
traction, are shown in Fig. 4(a); the fitted peaks are dis¬ 
played in panel (b). The peak splitting qj = {dh, ±Sk, 0) 
corresponds to 6h = 0.048(2) and 6k = 0.081(1), where 
the number in parentheses is the uncertainty of the last 
digit of the preceding number. [For comparison, fitting 
gaussian peaks while assuming a single domain yields po¬ 
sitions Sh — 0.059(2) and 5k = 0.074(2).] 

The rotation of the peaks away from 5h = 5k (that is, 
a rotation of the modulation away from the Cu-0 bond 



FIG. 4. (color online) Elastic intensities around the antifer¬ 
romagnetic wave vector at 1.5 K, (a) measured (with back¬ 
ground subtracted) and (b) fit, as described in text. 
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direction) is qualitatively consistent with the peak ro¬ 
tation seen previously in La 2 -xSra;Cu 04 with x = 0.12 
[34, 35] and in La 2 Cu 04_|_5 [15], although the magnitude 
of the rotation in the present case, 14(1)°, is considerably 
larger. [For the single-domain fit, the rotation would be 
7(2)°.] The intensity contribution from the (0,1,0) do¬ 
main is < 60% of that from the (1,0,0), which is roughly 
consistent with the intensity anisotropy observed in the 
spin-glass regime for LSCO with x = 0.05 [36]. It sug¬ 
gests a tendency for spins to align along h (but with disor¬ 
der), which is also indicated by the anisotropic magnetic 
susceptibility measured on detwinned crystals oix = 0.03 
[37]. In our fitting, the gaussian widths were allowed to 
be anisotropic along the h and k directions, resulting in 
ah = 0.045(3) rlu and ak = 0.028(2) rlu. 

The temperature dependence of the elastic peaks was 
probed in two ways. With Ef = 3.7 meV, a scan 
was measured along Q = (1.07, A:, 0) at each temper¬ 
ature. The intensity was fit with a pair of gaussians 
plus constant background. Both the integrated intensi¬ 
ties and the peak amplitudes are plotted in Fig. 5. With 
Ef = 5 meV, the intensity was measured at a single po¬ 
sition, (1.07, —0.07,0), as a function of temperature, and 
the measurements have been averaged over bins of 5-K 
width to improve statistics; background was obtained by 
assuming that no signal was present above 25 K. The 
resulting background-subtracted amplitude is also pre¬ 
sented in Fig. 5. 

While the statistical error bars are substantial, there 
appears to be a clearly distinguishable difference in the 
temperature dependence for the two sets of measure¬ 
ments. The onset temperature is correlated with the 
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FIG. 5. (color online) Temperature dependence of the inten¬ 
sity of elastic magnetic scattering, as described in the text. 
Filled symbols represent results obtained with Ef = 3.7 meV; 
(blue) squares for integrated intensity, (black) circles for peak 
amplitude; results normalized to 1 at the lowest tempera¬ 
ture. Open circles (red) indicate peak amplitude obtained 
with Ef = 5 meV. Error bars here and elsewhere correspond 
to la determined from counting statistics. 
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energy resolution for the measurement. The energy res¬ 
olution is slightly coarser ior Ef =5 meV, and the onset 
of intensity appears to be close to For the higher- 
resolution measurements with Ef = 3.7 meV, the main 
growth in intensity occurs below ^ 10 K. Going to yet 
higher energy/time resolution, muon spin rotation stud¬ 
ies indicate the onset of a finite hyperfine field at ~ 4 K 
[38, 39]. The resolution-dependence of the onset tem¬ 
perature suggests the gradual freezing of slowly fluctuat¬ 
ing spin correlations, as found in the spin-glass phase of 
La 2 _a;Sr 2 ,Cu 04 for x < 0.05 [40, 41]. 


2. Inelastic 


For low-energy spin fluctuations, the intensity has an 
incommensurate pattern similar to that of the elastic sig¬ 
nal in Fig. 4(a), though without the peak rotation and 
with somewhat greater peak widths, as originally demon¬ 
strated by Cheong et al. [42]. Here we focus on charac¬ 
terizing the frequency dependency of the magnetic sig¬ 
nal at base temperature (< 5 K ^ Tc) for the range 
0.2 meV < ?ia; < 15 meV. Figure 6 shows representative 
scans (from SPINS) and cuts (from HYSPEC data) for 
constant energies as a function of Q, as well as images of 
the intensity as a function oi E = Tiuj and Q. The lines 
through the data points in Fig. 6(e)-(f) indicate fitted 
gaussian peaks used to evaluate the integrated intensity. 
For each measurement condition, a distinct range of exci¬ 
tation energies is covered; in each case, there is substan¬ 
tial intensity at the lowest resolvable energies. For the 
gaussian fits, the peaks are symmetric in fc, with peak 
center and gaussian width allowed to vary independently 
for each set of data, but constrained to be independent of 
excitation energy. All of the results are consistent with 
a peak center of 5k = 0.060(4) rlu and gaussian width 
CTfc = 0.041(5) rlu. The magnitude of incommensurabil¬ 
ity for the inelastic scattering is about 10% smaller than 
that of the elastic signal; this is qualitatively consistent 
with the results for LBCO with x = 0.095 [9]. 

The integrated intensities obtained from the fits to the 
spin-fluctuation scattering shown in Fig. 6 are presented 
in Fig. 7. Here, each data set has been normalized to a 
measurement of the elastic incoherent scattering, which 
is dominated by the incoherent nuclear scattering from 
the sample; the spread in the results is consistent with 
the uncertainty in the normalization. We observe a peak 
at A « 0, with the intensity falling off to a constant level 
for A > 5 meV. The line through the data corresponds 
to a constant plus a lorentzian centered at if = 0 with 
a half-width of 1.3(2) meV. The absence of a spin gap 
in the superconducting state is consistent with a recent 
study of La 2 -a;Baa;Cu 04 with x = 0.095 [9]. 
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FIG. 6. (color online) Image plots of intensity vs. E and Q for (a) SPINS with Ef = 3.7 meV, (b) SPINS with Ef = 5 meV, 
(c) HYSPEC with Ei = 8 meV, (d) HYSPEC with Ei = 27 meV. Representative scans for (e) Ef = 3.7 meV, (f) Ef = 5 meV. 
Cuts integrated over 0.9 < h < 1.1 and —0.25 < I < 0.25 for (g) Ei = 8 meV, (h) Ei = 27 meV. Lines in (e)-(h) represent fits 
of a pair of symmetric gaussian peaks on a Q-independent background. Scans are vertically offset for clarity. In (c), there is a 
spurious feature (of undetermined origin) at E « 1.7 meV. 


C. Lattice response 

1. Phonons 

Soft phonons corresponding to octahedral tilts that 
break the equivalence of in-plane O sites were recently 
identified and studied in La 2 _a;Ba 2 ;Cu 04 with x = 
0.125 [23]; ordering of these displacements leads to a 
reduced lattice symmetry that can pin charge stripes 
[22, 43]. Given the evidence for charge stripe order in 
La 2 _a;Sr 2 ,Cu 04 [17, 35, 44, 45], we decided to test for 
such soft tilts. We did this by measuring about (3, —3,0), 
the same position studied in La 2 -a;Baa;Cu 04 and a wave 
vector at which elastic scattering is forbidden in the LTO 
phase. 



FIG. 7. (color online) Compilation of the Q-integrated mag¬ 
netic intensity (normalized as discussed in the text) as a func¬ 
tion of excitation energy. The line through the data is a fit to 
a constant plus a lorentzian centered at E = 0; the half-width 
at half-maximum of the lorentzian is 1.3(2) meV. 


Representative results, obtained aXT = 150 K, are pre¬ 
sented in Fig. 8(a)-(c). Panel (a) shows inelastic signal, 
integrated between 3 and 7 meV, and plotted vs. wave 
vector. The elliptical features at (2,—4,0), (2,—2,0), 
and (4,—2,0) are acoustic phonons about fundamental 
Bragg peaks, while the intensity at (3, —3, 0) corresponds 
to soft tilt modes. Panels (b) and (c) show the disper¬ 
sion through (3,—3,0) along longitudinal and transverse 
directions, respectively. One can see a continuous con¬ 
nection to transverse acoustic phonons in (c). For the 
longitudinal fluctuations in (b), the intensity appears to 
be interrupted by an anti-crossing with an unseen mode 
at ^ 12 meV that also impacts the longitudinal acoustic 
mode dispersing from (2,—2,0); nevertheless, the upper 
part of the the dispersion can be seen between 15 and 20 
meV. If the signal at (3, —3,0) were a soft mode, we would 
expect to see an energy gap; however, there is no sign of 
one. To the extent a gap might be present, it would have 
to be highly over-damped, which does not seem to agree 
with the dispersion clearly seen in Fig. 8(c). 

We repeated these measurements at several tempera¬ 
tures. To compare the temperature dependence of these 
data, we integrated the signal centered on (3, —3, 0) from 
3 to 7 meV; the results, after background subtraction, are 
indicated by red squares in Fig. 8(d). For comparison, 
the blue circles indicate a similar integral over acous¬ 
tic phonons about the (2, —2, 0) Bragg peak. The in¬ 
tensity at the “forbidden” position shows little change 
with temperature between 150 and 300 K; this is sim¬ 
ilar to the constant intensity found for the same mode 
within the LTO phase of La 2 _a;Ba 2 ,Cu 04 with x = 0.125 
[23]. On the other hand, there is some decrease in in¬ 
tensity on cooling to 80 and 5 K. The main observation 
here is that the two temperature dependencies are dif¬ 
ferent: that of the lattice phonon at the fundamental 
Bragg reflection is consistent with the detailed balance 
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FIG. 8. (color online) Measurements of phonons in La 2 - 2 :Sri:Cu 04 with x = 0.07 at T = 150 K. (a) Constant-energy slice 
{hoj = 5 ± 2 meV with L integrated from —0.25 to 0.25), showing strong acoustic phonons about the (2, —2, 0), (2, —4, 0), and 
(4, —2, 0) Bragg points, and soft tilt fluctuations at (3,—3,0). Dashed lines indicate the directions and widths of the slices 
through (3, —3, 0) shown in (b) (longitudinal direction) and (c) (transverse direction), (d) Intensity of the phonon signal at 
(3, —3, 0) integrated from 3 to 7 meV [as in (a)] and multiplied by 2.5 (red squares) and at (2,—2,0) (blue circles), corrected 
for background measured in between these positions. 


and therefore a T-independent imaginary part of the dy¬ 
namical susceptibility, x", while the tilt mode indicates 
temperature-dependent x" ^ with the characteristic tem¬ 
perature of ^ 150 K. 



FIG. 9. (color online) Elastic scattering (hcj = 0 ± 0.5 meV) 
about the (3, —3, 0) position along the (a) transverse and (b) 
longitudinal directions. In each case, the red curve is a gaus- 
sian fit to the (3, —3,0) peak. The extra peaks in (b) are 
powder diffraction peaks from the A1 sample holder. 


2. Superlattice peak 

Given the strong LTT fluctuations, we have also looked 
at the elastic scattering at the (3, —3,0) superlattice 
position. In-plane cuts along transverse and longitudi¬ 
nal directions are shown in Fig. 9; the out-of-plane cut 
is shown in Fig. 10. We find a peak that appears to 
have a resolution-limited width in all three directions. 
The temperature-dependence of the integrated intensity 
(evaluated from the transverse cuts) is shown in Fig. 11. 
It decays only gradually with temperature, and heads 
to zero close to the estimated temperature for the tran¬ 
sition to the high-temperature tetragonal (HTT) phase 
[46]. Given that the same superlattice peak was found to 
go to zero in the LTO phase of LBGO x = 0.125 under 
similar measurement conditions [23], it seems unlikely 
that the present peak could have a spurious cause, such 
as double scattering. Hence, it appears that the symme¬ 
try of the crystal structure is lower than that suggested 
by the analysis of neutron powder diffraction data [26]. 
As the average structure has a significant orthorhombic 
strain, the likely structure is the low-temperature less 
orthorhombic (LTLO; Pccn) phase. 

It is of interest to compare the intensity of the su¬ 
perlattice (3, —3, 0) peak to a fundamental Bragg peak 
such as (2,—2,0). At 5 K, we find an intensity ratio of 
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FIG. 10. Elastic scattering along Q = {3,—3,1). Based on a 
comparison with a cut along (2, —2, 1), we conclude that the 
peak shape and width are due to the vertical focusing and the 
sample geometry. 


7 X 10“^; this is an upper limit, as the (2, —2,0) intensity 
almost certainly saturated the detector, as well as suffer¬ 
ing secondary extinction from our large crystals. Such an 
intensity ratio is definitely outside of the dynamic range 
of a powder-diffraction measurement. We can also com¬ 
pare it with the ratio of low-energy phonon intensities, 
as evaluated in Fig. 8(d); the acoustic phonon intensity 
at wave vectors close to the Bragg peak should be pro¬ 
portional to the elastic peak intensity. From the figure 
we find an intensity ratio at 5 K of 0.4, three orders of 
magnitude greater than the elastic intensity ratio. This 
indicates that there is substantial entropy remaining in 
these tilt fluctuations at low temperature. 



FIG. 11. Integrated intensity of the (3,—3,0) superlattice 
peak as a function of temperature. The vertical dashed line 
indicates the transition to the HTT phase based on an inter¬ 
polation formula [46]. 


IV. DISCUSSION 
A. Glassy spin-stripe order 

The elastic magnetic intensity data in Fig. 5 suggest 
an onset temperature in the range of 10-20 K. A consid¬ 
erably lower transition temperature is indicated by nu¬ 
clear magnetic resonance (NMR) and muon spin rotation 
(/rSR) techniques, which are sensitive to much smaller 
time scales. ^^®La NMR and nuclear quadrupole reso¬ 
nance (NQR) studies of LSCO indicate a spin freezing 
temperature of 5 K for x = 0.06 [47] and a transition be¬ 
low 4.2 K for X = 0.07 [48]. This frequency dependence 
of the freezing transition indicates glassy behavior. In 
fact, the transition temperature increases smoothly on 
reducing the doping into the non-superconducting regime 
0.02 X ^ 0.055 [39, 49], where spin-glass behavior was 
originally identified in measurements of the bulk magne¬ 
tization [50, 51]. 

Even before the glassy response was experimentally de¬ 
tected, it was proposed that the doping of holes into anti¬ 
ferromagnetic Cu 02 layers would lead to a “cluster” spin 
glass due to frustrated electronic phase separation [52]. 
Indeed, /xSR measurements at T < 1 K have identified 
an average local hyperfine field of finite magnitude for 
0 < X < 0.09 [39], especially including x = 0.07 [38], 
consistent with local regions of antiferromagnetic order. 
The finite elastic peak width that we have observed in¬ 
dicates high sensitivity to quenched disorder; note that 
our doping level of 7% is far below the percolation limit 
for magnetic dilution [53]. In the context of coupled spin 
and charge stripe order [43, 54], it was proposed that the 
charge stripes would be the component that would couple 
strongly to electronic disorder, resulting in a stripe glass 
[55]. The generally strong sensitivity of charge-ordered 
phases to disorder in quasi-2D layered materials, such as 
cuprates, is also expected on the basis of the low-energy 
mapping of the interactions responsible for the ordering 
onto an effective random-field Ising model [56]. Neu¬ 
tron scattering measurements eventually demonstrated 
the presence of incommensurate elastic magnetic scat¬ 
tering throughout the spin-glass regime [12], consistent 
with a stripe glass phase. (An incommensurate spin- 
spiral phase has also been proposed [57]; however, the 
nearly uniaxial spin anisotropy that we observe is not 
consistent with a spin-spiral state.) 

While the absence of long-range stripe order in the 
presence of disorder is inevitable in two dimensions, a 
new possible glassy state, the spin-density-wave (SDW) 
glass, has recently been proposed [58]. In the SDW glass, 
the spins maintain a preferred direction over long dis¬ 
tances (spin nematic order). It is the proliferation of 
single dislocations of the charge stripes [59] that would 
normally destroy such order. It is presumed that, for 
sufficiently weak disorder, a regime may exist in which 
double dislocations proliferate but single dislocations are 
suppressed. 

Could our sample exhibit SDW-glass behavior? The 






fact that we observed the rotation of the incommensu¬ 
rate magnetic peaks away from the Cu-0 directions in 
the presence of structural twin domains is a consequence 
of the tendency for the spins to align along [010], as dis¬ 
cussed in Sec. IIIBl. Of course, this is essentially the 
same orientation of the spins found in antiferromagnetic 
La 2 Cu 04 , where the spin orientation is determined by de¬ 
viations from pure Heisenberg spin coupling due to weak 
exchange anisotropy resulting from spin-orbit coupling 
effects [60]. Hence, while we do see a tendency towards 
spin-nematic order (and away from a fully random distri¬ 
bution of spin orientations), it is not clear whether this 
is direct evidence for the SDW-glass state. 


B. Pair-density-wave glass 

We have observed quasi-static incommensurate anti¬ 
ferromagnetism, with gapless magnetic excitations, at 
temperatures well below the superconducting transi¬ 
tion. These results are similar to those found in 
Lai,go 5 Bao.o 95 Cu 04 , where weak charge and spin stripe 
order [61] and gapless spin excitations [9] coexist with 
superconductivity. Gapless incommensurate spin fluctu¬ 
ations have been seen previously in LSCO with x ^ 0.13 
[19, 20, 62, 63], although this behavior has not always 
been recognized as intrinsic for measurements without 
an applied magnetic field [63]. 

As noted elsewhere [9], the coexistence of quasi-static 
spin stripes with superconductivity is inconsistent with 
expectations for a uniform d-wave superconductor [3, 4] 
but is compatible with a spatially-modulated supercon¬ 
ducting order parameter, as with the proposed pair- 
density-wave state (PDW) [2, 10, 64, 65]. The enhance¬ 
ment of spin-stripe order in underdoped La 2 -a;Sr 2 ;Cu 04 
by application of a c-axis magnetic field has been ob¬ 
served previously [66]; such a field also leads to decou¬ 
pling of the superconducting layers [67, 68], the effect 
that originally motivated the concept of the PDW super¬ 
conductor [10, 64]. 

A concern with the PDW state has been that it is very 
sensitive to disorder [11]. A distribution of single disloca¬ 
tions can destroy the long-range phase order. It happens 
that the SDW glass theory also applies to a PDW glass 
[58]; in the latter case, it is proposed that superconduct¬ 
ing order may survive in the form of a 4e nematic phase 
[69]. (The transformation from PDW order to 4e nematic 
superconductor due to thermal fluctuations is discussed 
in [70].) 

Another potential challenge to obtaining supercon¬ 
ducting order from PDW correlations is the frustration 
of the interlayer Josephson coupling, as occurs in he case 
of LBCO with X = 0.125 [64, 71]. In the present case, 
the crystal structure is orthorhombic and the stripes are 
rotated from the Cu-0 bond direction, so that a finite 
interlayer Josephson coupling is allowed (and observed 
[72]). Hence, the presence of a PDW superconducting 
glass state in our LSCO sample is quite plausible. 


C. Inequivalent oxygen sites and stripe pinning 

The onset of charge and spin stripe order in LBCO is 
limited by a transition to the LTT (or LTLO) structure 
that results in inequivalent oxygen sites within the Cu 02 
planes [16], as illustrated in Fig. 1. Hence, it has been 
surprising that in LSCO, with a presumed average LTO 
structure and corresponding equivalency of all in-plane 
O sites [26], the NMR evidence for slow spin and charge 
fluctuations begins at temperatures even higher than in 
LBCO [17, 73, 74]. While our evidence for a lowered 
LTLO symmetry resolves this puzzle, it is of interest to 
consider relevant experimental results from the past. 

In an early study of LSCO with x = 0.14, Nohara et al. 
[75] found a softening of the transverse elastic constant 
(cii — ci 2)/2 below 50 K. They noted that the associated 
strain on the lattice would lead to inequivalent oxygen 
sites, as in the LTLO structure. Later tests of struc¬ 
ture tended to focus on a; = 0.12. Transmission electron 
microscopy (TEM) studies demonstrated the presence of 
(llO)-type superlattice peaks, consistent with either the 
LTLO phase or LTT phase [76-78]. Imaging these peaks 
revealed that the lower-symmetry phase is present at the 
domain boundaries between twins of the LTO phase. A 
synchrotron x-ray powder diffraction study by Mooden- 
baugh et al. [79] found evidence for a small fraction 
(~ 10-20%) of LTT phase within the LTO phase of LSCO 
x=0.12 for T < 100 K [79]. For this same composition, an 
NMR study reported a change in symmetry below 50 K, 
consistent with LTT-like tilts [80]. 

In the present study, as the neutron diffraction mea¬ 
surement averages over the sample volume, we have to al¬ 
low the possibility that the weak (330) superlattice peak 
could be coming from a small fraction of the sample, 
such as twin boundaries. At the same time, the strength 
of the octahedral-tilt phonons indicates a response from 
the entire volume. That a small fraction of these fluctua¬ 
tions might condense throughout the volume should not 
be surprising. In fact, there might be a small tilt order¬ 
ing in the bulk and a much larger ordering in the twin 
domains. A local-probe measurement, such as TEM, will 
be necessary to provide a definitive answer. 

Previous studies of soft phonons in LSCO [81, 82] have 
focused on the octahedral tilts about plaquette diagonals. 
In the high-temperature tetragonal (HTT) phase, there 
are two such modes that are orthogonal and equivalent. 
One of these effectively goes soft at the transition to the 
LTO phase [21], resulting in superlattice peaks of the 
type (032). The possible softening of the other mode can, 
in principle, be studied at positions such as (302); how¬ 
ever, the inevitable presence of twin domains means that 
the (302) of one domain is nearly superimposed on the 
(032) of its twin, complicating the measurements. On the 
other hand, a position such as (330) is a forbidden super¬ 
lattice reflection in the LTO phase, but is allowed in the 
LTT and LTLO phases [16]. As already mentioned, soft 
phonons at this position directly correspond to displace¬ 
ments that cause in-plane oxygens to be inequivalent. In 
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a recent study of LBCO with x = 0.125 [23], soft phonons 
were observed at (330) with an integrated intensity that 
was independent of temperature in the higher tempera¬ 
ture phases (HTT and LTO); on entering the LTT phase 
at 54 K, the intensity of the soft phonons was transferred 
to the superlattice peak at zero energy. 

For LBCO, theory finds, correctly, that the ground 
state has static LTT tilts, but that the LTO structure 
is favored at higher temperatures due to the free-energy 
gain from the higher entropy of that state [83, 84]. The 
LTT tilts still cost little energy, but contribute entropy as 
they are unable to order in the LTO phase. For our LSCO 
X = 0.07 sample, we appear to have an intermediate situ¬ 
ation: a small fraction of the LTT-like tilt character con¬ 
denses, but most of it is prevented from doing so by the 
orthorhombic symmetry, which likely becomes LTLO in¬ 
stead of LTO. The LTT-like tilt fluctuations have a quite 
noticeable impact on the thermal conductivity. In the 
case of LBCO, complete or nearly complete ordering of 
the tilts results in a substantial increase in the thermal 
conductivity below the transition temperature [31]. In 
contrast, thermal conductivity measurements on LSCO 
with X > 0.1 show no significant recovery at low temper¬ 
ature [85, 86], consistent with the continued presence of 
dissipative fluctuations. 

The impact of the LTT-like phonons should have only 
an indirect impact on the spin fluctuations. Condens¬ 
ing some of these phonons to induce the LTLO phase 
allows charge stripes to be pinned, which in turn allows 
the spin stripes to become quasistatic. In contrast, a dif¬ 
ferent phonon mode has recently been shown to couple 
directly to the spin fluctuations [87, 88]; it involves in¬ 
plane displacements of O atoms in directions transverse 
to the Cu-0 bonds [87]. 

It is interesting to note that Baledent et al. [89] studied 
a weak (110) peak, equivalent to our (330), that appeared 
below 120 K in an LSCO crystal with x = 0.085 and 
Tc = 22 K. They observed this feature with polarized 
beam in the spin-flip channel, and hence interpreted it 
as a magnetic diffraction peak. Given the present results 
regarding structural correlations, it may be worthwhile 
to test the degree of magnetic character at (110). 


D. Rotation of incommensurate peaks 

The rotation of the IC magnetic peaks away from the 
Cu-0 bond direction is consistent with previous measure¬ 
ments on LSCO X = 0.12 [34, 35] and La 2 Cu 04+5 [15]; 
it is also consistent with the orthorhombic symmetry, as 
discussed in [90]. The magnitude of the rotation found 
here for x = 0.07, ~ 14°, is significantly larger than that 
reported for x = 0.12 (^ 3°). For x ~ 0.055, the rota¬ 
tion reaches 45° [12], with the spin modulation occurring 
uniquely along the b axis in the spin-glass regime [36]. 

The rotation of the modulation direction at x ~ 0.055 
corresponds with the superconductor to insulator transi¬ 
tion. A careful study of resistivity in an LSCO thin film 


with carrier concentration varied electrostatically pro¬ 
vides evidence that the transition involves a localization 
of pairs [91]. A torque magnetometry study found ev¬ 
idence that superconducting fluctuations survive on the 
insulating side of the transition [92]. On the larger x side, 
measurements of hysteresis in magnetoresistance indicate 
charge glass character of the insulator state, with resid¬ 
ual effects observable in a superconducting x = 0.06 thin 
film [93]. Even in zero field, there is upward curvature 
of Pab{T) at T < 100 K for X < 0.13 that grows as x 
decreases toward 0.06 [94]. 

Hence, it appears that the rotation of the IC peaks in 
our X = 0.07 sample may involve a mixing of the charac¬ 
ters, superconducting and insulating, respectively associ¬ 
ated with the bond-parallel and plaquette-diagonal mod¬ 
ulations. The insulating character need only involve lo¬ 
calization, rather than the breaking, of pairs. As doping 
increases, the proportion of diagonal/insulating charac¬ 
ter decreases, as indicated by the reduced rotation angle 
found for X = 0.12 [34]. 


E. Connection with nodal gap? 

The presence of a spatially-modulated superconduct¬ 
ing state should have an impact on the observed super¬ 
conducting gap. When the superconductivity is uniform, 
the gap is observed to have dx'i_y 2 symmetry, with gap¬ 
less nodes along the (1,1) and (1,-1) directions [95]. 
Interestingly, anomalous behavior has been detected in 
LSCO with X ^ 0.1. In a recent angle-resolved photoe¬ 
mission spectroscopic (ARPES) study of LSCO x = 0.08 
[96], the appearance of a gap at the nodal points was in¬ 
ferred, starting at temperatures as high as ~ 80 K and 
growing on cooling. Given that the resistivity on such 
a sample does not show either strong insulating or su¬ 
perconducting behavior at such high temperatures, the 
nature of the gap is somewhat ambiguous; nevertheless, 
there is a clear change in the spectral function at the 
nodal wave vector as a function of temperature. 

To put this in perspective, we can compare with the 
case of LBCO with x = 1/8, where the presence of PDW 
order has been inferred [64, 71]. Theoretically, the PDW 
superconductor is predicted to have a large gap in the 
antinodal regions, but to have no gap along an arc of 
states centered on the nodal direction [97, 98]. Experi¬ 
mentally, a d-wave like gap is seen in ARPES measure¬ 
ments about the nodal regime [99, 100]; the deviation 
from the PDW prediction could be due to the presence 
of spin stripe order, which was not included in the theory 
for the gap [2]. 

Further evidence that the gap in underdoped LSCO 
is different from that in LBCO has been provided by 
studies of thermal conductivity. At low temperature, the 
thermal conductivity, k, is dominated by the contribu¬ 
tion from electrons, and in the presence of a d-wave gap, 
only quasiparticles near the nodes can contribute. It has 
been shown that the ratio of the thermal conductivity 
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to temperature, n/T, in the limit T —)■ 0, has a universal 
form that depends only on the Fermi velocity and the gap 
velocity near the node [101]. Low-temperature measure¬ 
ments of K have provided measures of the superconduct¬ 
ing gap in a variety of cuprates that are consistent with 
results from ARPES [102]; however, anomalous behavior 
is observed in LSCO with a; < 0.14 [102-104]. Effectively, 
the thermal conductivity is smaller than one would ex¬ 
pect from well-defined quasiparticles, and the magnitude 
decreases in a magnetic field, contrary to conventional 
behavior. 

One way to explain the results for underdoped LSCO 
is to introduce a small gap at the nodal point [105]; how¬ 
ever, to explain the response to a magnetic field, one must 
then assume that the gap grows with field. Alternatively, 
it has been proposed that the field induces local antifer¬ 
romagnetic order within the cores of the field-induced 
vortices, thus suppressing the quasiparticles essential for 
heat transport at low temperature [106]. The case of 
LBCO with X = 1/8 is an example where local antifer¬ 
romagnetic order, in the form of spin stripes, does not 
result in a nodal gap [16, 100]. A difference in the case of 
LSCO is the rotation of the stripe orientation. The mix¬ 
ing in of diagonal stripe character might be responsible 
for the apparent nodal gap [107]. 


V. SUMMARY 

We have presented a neutron scattering study of LSCO 
X = 0.07, an underdoped cuprate with bulk superconduc¬ 
tivity and Tc = 20 K. We have observed incommensurate 
spin fluctuations with a gapless spectrum that coexists 
at low temperature with the superconductivity. These 


results are consistent with a periodic modulation of the 
superconducting state intertwined with the spin order 
[2], and they demonstrate that the related behavior pre¬ 
viously observed in LBCO [8, 9, 71] is not unique. We 
have also discovered evidence of a lowering of symme¬ 
try below that of the presumed LTO phase, likely to 
the LTLO; in addition, considerable entropy remains in 
LTT-like phonons at low temeperature. These observa¬ 
tions provide a connection to the stripe pinning found in 
the LTT phase of LBCO [16]. The weakly-ordered spin 
stripes are rotated away from the Cu-0 bond directions, 
consistent with previous studies of LSCO with x = 0.12 
[34, 35], and suggesting pinning of charge stripes in the 
proposed LTLO phase. We have pointed out that the 
mixing in of some diagonal stripe character may help to 
explain a recent ARPES observation of a nodal gap that 
appears even above Tc [96]. 
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